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Abstract 

'Trt ' We present here an overview of the Fourier Transform Scanning Tunneling spectroscopy technique 

rS^ . (FT-STS). This technique allows one to probe the electronic properties of a two-dimensional system 

c/3 ! by analyzing the standing waves formed in the vicinity of defects. We review both the experimental 

and theoretical aspects of this approach, basing our analysis on some of our previous results, as well 
as on other results described in the literature. We explain how the topology of the constant energy 
C^ ■ maps can be deduced from the FT of dl/dV map images which exhibit standing waves patterns. We 

show that not only the position of the features observed in the FT maps, but also their shape can 

i-A . be explained using different theoretical models of different levels of approximation. Thus, starting 

C I with the classical and well known expression of the Lindhard susceptibility which describes the 

screening of electron in a free electron gas, we show that from the momentum dependence of the 
susceptibility we can deduce the topology of the constant energy maps in a joint density of states 
approximation (JDOS). We describe how some of the specific features predicted by the JDOS are 

^T . (or are not) observed experimentally in the FT maps. The role of the phase factors which are 

'!;;j- I neglected in the rough JDOS approximation is described using the stationary phase conditions. 

^;5 ! We present also the technique of the T-matrix approximation, which takes into account accurately 

these phase factors. This technique has been successfully applied to normal metals, as well as to 

;<— ^ . systems with more complicated constant energy contours. We present results recently obtained 

on graphene systems which demonstrate the power of this technique, and the usefulness of local 
measurements for determining the band structure, the map of the Fermi energy and the constant- 
energy maps. 
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INTRODUCTION 

One of the most remarkable feats achieved with an STM, beside the possibihty to visuahze 
material surfaces with atomic resolution, was the possibility to image the standing waves 
associated to the interference of quasi-free electron wavef unctions. This has been achieved 
for the first time on a copper surface, for electrons confined in a circular resonator created 
with iron ad-atoms, structure well known as a "quantum coral" \Mi^- This observation 
has provided direct evidence that electrons are associated to waves, thus demonstrating the 
wave-particle duality which is one of the fundamental concepts of quantum mechanics. 

The "standing waves" arising in the presence of surface inhomogeneities are also known 
as Friedel oscillations [5] , where the term of Friedel oscillations has been first introduced to 
describe the asymptotic dependence of the perturbed density of states of a free electron gas 
in the presence of disorder. Their observation allows the illustration of some very important 
concepts of condensed matter physics. Thus, the analysis of Friedel oscillations provides a 
direct observation of screening and of electron-electron interaction. Moreover, these oscil- 
lations lie at the foundation of the description of the indirect coupling between magnetic 
moments via the conduction electrons in a metal with the famous Ruderman-Kittel-Kasuya- 
Yosida (RKKY) interaction potential p-lsl , as well as of the long range adsorbate interaction 
mediated by a two-dimensional electron gas [9| . 

Following their first description by Friedel, the possibility to use the Friedel oscillations to 
probe the electronic structure of materials was considered by many others, notably in the case 
of transitional metals [lO|. They were mentioned in relation to magnetic impurities in bulk 
(3D) materials, for which an important damping factor of the amplitude of the oscillations 
is observed: the oscillations fall off with the distance r as 1/r" where a is the dimensionality 
of the considered electron gas. Furthermore, the observation of the Friedel oscillations was 
for the first time done indirectly by the observation of the coupling between two magnetic 



layers with a non-magnetic spacer (host metal) [ll|| for which the theoretical description [12 
has revealed the Fermi surface of the spacer. The development of the scanning tunneling 
microscopy (STM) has offered the possibility to study the standing waves in the local density 
of states (LDOS) which are in fact energy resolved Friedel oscillations for 2D or ID electron 
gases which provide longer coherence lengths than those observed in 3D. The local density 



of surface states has been first obtained by Hasegawa and Avouris [13[ on confined Shockley 
states of Au(lll) surfaces at room temperature, and by Crommie et al. 1^ on Cu(lll) at 
4K. 

Subsequently, the dispersion relation E{k) of surface-state electrons of Ag(lll) and 
Cu(lll) permitted the estimation of the surface state inelastic lifetime [l5|. As large en- 
ergies are accessible both below and above the Fermi level, it was also possible to study 
the deviation of the free-electron-like parabolic dispersion, moreover Biirgi et al. [l6| have 
defined a method to directly image the potential landscape on Au(lll) by STM. 

For these measurements, the standing waves have been imaged not only at the Fermi 
level but also at different energies. This requires the spatial mapping of dl/dV at different 
applied voltages V, which allows one to focus on individual values of the energy specified by 
V, and avoids an integration over all the wavelengths corresponding to the energies between 
the considered energy level and the Fermi level. This was clearly demonstrated by Schneider 



et al. [17| on Ag(lll). In this study, standing waves are generated by step edges, and one 
focuses on a particular reciprocal lattice direction when analyzing the Fermi surface. 

In 1997 Sprunger et al. jl8| demonstrated the possibility to use STM to directly image 



the Fermi surface contour of a metal surface, by performing the power spectrum (Fourier 
transform (FT)) of a topographic image of a complex "electron sea" pattern generated by a 
random distribution of point-like surface defects. These studies were performed around the 
Fermi level (low bias voltage) and on simple isotropic Fermi surface of Au(lll) and Cu(lll) 
surfaces [19|. In this seminal paper Friedel oscillations appear in the Fourier transform as 
a circular feature centered around the Brillouin zone center with a radius corresponding to 
2k f. This corresponds to the simple case of an isotropic Fermi contour centered around the 
center of the Brillouin zone. As we will see in the next section, this is a direct evidence 
of the singularity in the Lindhard susceptibility function of the two dimensional electron at 
scattering momentum vector q = 2 k j. 

The dependence of the Lindhard susceptibility on q was subsequently elegantly demon- 
strated by the observation of standing waves on Be (1010) surfaces |20|. In this system 
the Fermi surface is no longer isotropic and standing waves can not be generated by all 
step edges (i.e. in all directions). As we will see later this corresponds to the fact that 
the wavevectors for the quasiparticles with a given energy can isotropically point in any 
direction for a circular Fermi contour centered around the F, but acquire a dependence on 
the direction and position in momentum space if the Fermi surface is no longer isotropic, or 
if it becomes spht. 

More recently the FT-STM technique (obtaining and analyzing the Fourier transform of 
STM images) has begun to be applied also to LDOS images {dl/dV maps) which combine 
STM imaging and a spectroscopic measurement of the LDOS as a function of energy and po- 
sition. This technique has been denoted FT-STS: "FT-Scanning Tunneling Spectroscopy", 
though in the literature the terms FT-STM and FT-STS are sometimes both used to describe 
FT-STS measurements. The FT-STS technique has been applied on high-Tc superconduc- 
tors, for a small range of energy around the Fermi level (see e.g. |2l|-|23l|). Another important 



observation made using the FT-STS was about the spin of the quasiparticles. As two waves 
with opposite spin directions cannot generate quasiparticle interferences (QPI), Pascual et 



al. [2J] demonstrated the ability to probe indirectly the orientation of the spin associated to 
a Fermi surface of a half magnetic material. We have tested the strength of this technique 
on a semi- metal ErSi2, for which some constant energy contours (CEC's) are split into 
several bands, and we have demonstrated as well the possibility to determine the whole 2D 
band structure in a wide range of energy 25|, |26| . We have shown that the power spectrum 



features can be easily explained on the basis of a joint density of states (JDOS) approach, 
by a simple geometrical construction formally established in [271 ]. 

Here we provide a background for the FT-STS technique by reviewing some previous 
theoretical and experimental results. Thus, besides the JDOS approach, we describe also 
how the phase factors can be taken into account in describing the FT-STS features. We 
underline that the stationary phase conditions lead to some selection rules in the scattering 
events, as considered in the case of the Kohn anomalies developed for example by Roth 



et al. [28| for the theoretical determination of RKKY interactions for non-spherical Fermi 
surfaces. A more accurate and complete theoretical T-matrix calculation is also presented. 
We show then how this technique applies to epitaxial graphene on SiC Furthermore, the 
determination of the band structure and of the Fermi velocity of graphene quasiparticles, 
the possibility to identify the position of the impurities, as well as the possibility to predict 
a large extension of the Van Hove singularity in epitaxial graphene with intercalated gold 
clusters are discussed. 



EXPERIMENTAL TECHNIQUE 



Our experiments were performed with a LT-STM from Omicron at 77 K at a base pressure 
in the 10~^^ mbar range. The dl/dV images were acquired using a lock-in amphfier and 
a modulation voltage of ±20mV. The graphene samples were prepared in UHV by the 
annealin g o f n-doped SiC(OOOl) at 900 K for several hours and subsequent annealing at 



1500 K |29H32|. This preparation method leads to the formation of a buffer graphene 
layer covalently bonded with the substrate and a monolayer graphene decoupled from the 
substrate [33|. The epitaxial graphene has an intrinsic n-type character and the Dirac point 



is at 0.4 eV below the Fermi level 134-36 



The deposition of gold on graphene was carried 
out at room temperature using a homemade Knudsen cell calibrated u sing a Quartz Crystal 
Microbalance. The sample was further annealed at 1000 K for 5 min |36| . 



FT-STS MEASUREMENTS ON AU(lll) AND ErSi2 AND THEIR INTERPRETA- 
TION BY THE JOINT-DENSITY-OF-STATES APPROXIMATION 

Background 



The underlying principle of the FT-STS technique stems from the screening of electrons 



around a localized impurity. As described in textbooks J37|, l38| , when a positive charge is 
held in a free electron gas, this charge will attract electrons, creating a surplus of charge 
which screens its electric field. It is common to solve the Poisson equation for this system 
in the momentum space. The response of the system can be described by the dielectric 

constant, which at a given momentum g is reduced to e{q) = 1 ^^^^, where q = \q\ and 

the x(Y) susceptibility is g -dependent in the reciprocal space. The susceptibility can be 
calculated using first order perturbation theory. In this approximation we consider that all 
eigenstates are "mixed" in the scattering process. Thus the Lindhard theory applies, and 
the susceptibility is given by [3j 



x{q) = Y1 



/(fc)-/(|fc + g1) 



(1) 



where f{k) is the Fermi-Dirac distribution function. The summation is performed over 
all possible k vectors providing that the state is occupied, as enforced by the Fermi-Dirac 
distribution. In the limit of vanishing q one obtains the Thomas-Fermi approximation, 
while for larger values of g, comparable to the value of the Fermi momentum at T = O-R', 
the susceptibility can be calculated explicitly, the result depending on the dimensionality of 
the electron gas. For a two-dimensional electron gas the susceptibility is given by: 



x'^{q)=n'^{Ep) 



1- 



q 



2kf 



e{q - 2kF) 



(2) 



where ri^^^Ep) is the two dimensional density of states at the Fermi level. The susceptibility 
is no longer analytic for q = 2k p. Here we have considered a nearly-free electron gas 
{E{k) = h'^k'^/2m*), for which the energy-resolved electron density oscillates with a wave- 
vector length q 



2't. 



Measurements of the FT-STS of the LDOS in Au(lll) and ErSi2 

A first experimental illustration of the FT-STS technique is described in figure [1] for 
Schockley states probed on Au(lll) surface. In a) the topographic STM image shows the 
herringbone reconstruction, and many defects (scatterers) are visible on the surface. The 
dl/dV map taken with a lock-in amplifier allows to measure directly the energy-resolved 
local density of states at several energies. We do not limit our analysis at the Fermi energy, 
but we analyze several different energies. When we perform the Fourier transform (here the 
power spectrum) of the different dl/dV images at increasing bias voltage we observe (in c)) 

a circle which increases in size. The radius of the circle corresponds to the | g |= 2 | k | as 
shown in the geometrical construction e). We report the value of g as a function of energy, 
which yields directly the dispersion curve, which is here parabolic, as expected in the case 
of a quasi "free" electron gas. This is a direct proof of the susceptibility singularity at 
momentum | g |= 2 | k |. 

This construction can be generalized for non circular CECs. In the right part of figure 
1, for comparison we show the FT-STS for a two dimensional electron gas of ErSi2- The 
system investigated here is the Erbium disilicide grown on Si(lll)-7x7. This system has 



been extensively studied theoretically and experimentally by photoemission [39|, |40|, and an 
atomic model of this structure has been determined by Gewinner and co-workers 39|, |41 . 
This structural model consists in a p(lxl) plane of Erbium atoms inserted between the 
silicon substrate and a buckled silicon layer J40|. This system leads to a perfect 2D metal- 
semiconductor interface which shows a semi-metallic character with a hole pocket centered 
around the center of the surface Brillouin zone F and six electron pockets around the M 
symmetry points at the Fermi level. As shown in figure [Hi) the expected dispersion is "hole- 
like" because it shows a reversed parabolic dispersion indicating a " negative" effective mass. 
Indeed, contrary to the case of Au(lll) surface, the size of the feature in the FT decreases 
when approaching the Fermi level. 

In f ) and g) we show typical topographic and dl/dV map images with point defects and 
standing waves pattern. Here the waves are non longer circular but have a hexagonal shape. 
In h) the power spectrum shows clearly a hexagon whose size decreases when approaching 
the Fermi level. The direct consequence of this CEC topology and of the fact that the 
susceptibility is q dependent, is the reinforcement of intensity along specific directions. 
Indeed along VM the FT is brighter than along TK. While at first glance it may seem 
that this could be an effect of the form of the matrix elements in the scattering process, it 
can however be explained entirely by the topology in the CECs. As schematized in figure 
[T]i), the number of q vectors corresponding to the scattering events at the susceptibility 
singularity can be symbolized by the area of the intersection between two CECs translated 
by g = 2k. For a circular contour this area is the same in all directions of the momentum 
space (as schematized in e) ), while for the hexagonal one this area is bigger in the VM 
direction (blue arrow) than in the TK direction (red arrow). 

The joint density of states approximation 

In a previous paper [271] we have formally established that the FT of a standing waves 
pattern image could be approximated by a joint density of states calculation (JDOS). In 
the presence of defects, the quasiparticles in the Bloch state k can be scattered into the 



Bloch state k , and new eigenstates can be constructed, in perturbation theory, as linear 
combinations of the degenerate unperturbed states k that belong to the constant energy 
contour E{k) = u. The leading term in the Fourier component of the LDOS at wave vectors: 
q = k' — k + G, where G is a reciprocal lattice vector, has an amplitude which takes the 
general form 



9{^,Q) 



^ fff{k,k',G)5iq~k + k'±G)£kd^k' (3) 



47r 

E{k)=E{k') 



where f{k,k',G) is a weighting factor that depends on scattering matrix elements, and 
thus on the overall distribution and nature of defects. 

The function g{u, q) gives the joint density of states with a main contribution for G = 

and replicas shifted by G . The quantity g{u, q) may be calculated by solving the Schrodinger 
equation for simple-defect geometries. Yet, if many defects of various symmetries are present, 
one may assume that f{k,k',G)isa fairly smooth, slowly varying function of k and k' . 
Thus, practically in the JDOS approximation, g{uj, q), which is related to the power spectrum 
of the LDOS, could be calculated by performing the self correlation function of the CEC at 
a given energy. In practice, the JDOS calculations consist in counting the number of pairs 
of k ,k' wave vectors, yielding a scattering wavevector g with the same length and direction. 
We fabricate an image where the pixels position is defined by the wavevector q, and the grey 
contrast is proportional to the number of the corresponding wavevectors pairs k, k'. This is 
a simple phenomenological approach of the gene ralization of the Lindhard susceptibility to 
Bloch waves states as developed by Blandin |42| . 

The figure [2] illustrates the fact that not only the structure, but also the shape of the 
features observed in the power spectrum, could be interpreted with a simple JDOS approach. 
In a) and b) are represented the CECs, the measured and theoretical band dispersion of the 
ErSi2 system respectively. Up to -250meV (rectangular red box in b)), the CEC is split 
into two types of contours, the hole-like band structure previously discussed in Figil] and, 
crossing the Fermi level from unoccupied states, six ellipsoidal electron pockets emerging 
around the M points. In the JDOS interpretation of the FTs, we should consider both 
the scattering processes joining hole-like and electron-like pockets (red arrow in a)), as well 
as scattering between the ellipsoidal electron-like pockets (blue arrow). Both should yield 
features centered around the M points. 

By changing the bias voltage in this energy range, the position of the feature does not 
change but only the shape. In c) we show the FT obtained at the Fermi level. One of the 
features arising in the FT is highlighted, and we can see that it has a "butterfly" shape. 
We have found that this particular shape is associated to the self correlation between two 
CECs, a central circle (the hole band) and an elliptic electron pocket, i.e. it is the result 
of the scattering events schematized by the red arrow in a). The JDOS calculation for 
electron-electron scattering is also showed in c). This clearly demonstrates that the leading 
scattering event joins hole-like and electron-like pockets. 

The dispersion relation as obtained from the FT-STS measurements is reported in b). 
This relation was deduced from the modification of the feature size with the bias voltage 
26( 1 . An excellent agreement between our experimental points and previous ARPES mea- 
surements 40(] is found. The agreement is also excellent with the calculated band structure 
from Rogero et al. 43| in the direction TM, but not in the direction TK. Obviously, the 



YSi2 DFT-LDA calculation reproduces the global shape of the measured hole band but the 



higher-energy excitations are shifted by values as large as 250 meV. Similar shifts for the 
predicted surface-state energy positions around the K point have been mentioned before for 



YSi2 and GdSi2 and the necessity to improve LDA self-energy term has been noted J43 



The case of similar dispersing CECs for high Tc superconductors has been studied in 



detail (see e.g. J22|, |23|), and more recently, for Bi2Sr2CaCu20s+s{Bi2212), the relation- 
ship between the ARPES measurements and the quasiparticule scattering interpretation of 
FT-STS measurement has been done using the concept of JDOS |44J. |45| . A more thorough 



theoretical approach to interpret these results, based on a Green function formalism, has 



been proposed by Wang and Lee |46|. Thus one can link the two measurement techniques 
yielding the real space single particle spectra for the STM and LDOS measurements, and 
the momentum space single particle spectral function for the ARPES measurements J45|. 
To summarize, as the first order stationary perturbation theory is invoked in the JDOS 
approach, and as large signal-to-noise ratio is expected with this technique it is crucial to 
apply it to physical systems in which the scattering processes mix sufficiently the unper- 
turbed wave functions in k-space and in energy. It has been discussed in detail by Capriotti 



et al. [47J. b v Kodra and Atkinson [48| , and for the case of high Tc superconducting mate- 



rials [211, |22| . We note that this technique has been successfully applied even more recently 
for G'aA^(OOOl) |49|. 

Beyond the JDOS approximation - the stationary-phase approximation 

Besides the electron-hole scattering processes, one should also take into account the hole- 
hole and electron-electron scattering events. At first glance, the predominance of electron- 
hole scattering events seems to be the effect of hidden matrix elements included in the 
function f{k,k',G) in the expression of the JDOS (equation E]). However we can show 
that a simple interpretation can already be provided using the background theory of the 
susceptibility. 

Up to now we have considered the steady states and neglected the phase factors in the 
rough JDOS approximation. In order to take into account these phase factors, the method 
of the stationary phase approximation is usually applied. It consists in noting that in the 
integral of the JDOS formula, oscillatory term with rapidly-varying phase will cancel, while 
the terms with the same phase should be added together. This leads to a supplementary 
criterion in the geometrical construction of the JDOS, similar to the case of the Kohn 



anomaly. As illustrated in figure Oa) and along the lines of Ref. |42[, we consider that the 



leading scattering process take place for values of the (f momenta joining points of the CECs 



for which Vfc(-E(/c)), 'Vk'{E{k')) (or alternatively the tangents to the CECs) are parallel and 
point in the same direction. 

In b) we consider a CEC similar to the Fermi surface of the Erbium disilicide with the 
central circle and the six ellipses. The full calculation of the JDOS is given in c) where 
we can easily recognize the features associated to the hole-electron and electron-electron 
scattering processes. By applying however the above-described criterion in our calculation 
we have obtained the JDOS given in d). Only the "butterfiy" features are preserved in 
the resulting spectra. This shows the importance of the phase factor, as well as that the 
matrix elements have no real effect on the FT features. This experimental result also nicely 
demonstrates the "stationary phase" mathematical theorem. 



FT-STS MEASUREMENTS ON GRAPHENE AND THE T-MATRIX APPROXI- 
MATION 

T-matrix approximation 

The T-matrix approximation has been successfully used to calculate the effects of disorder 
on the spectral properties of an electronic system. The basic theory is described for example 
in J50|, l5l|, and for a few example of how this is applied we mention Refs. 52|-|70| . The 



basic principle of this technique lies in an infinite perturbative summation of the diagrams 
resulting from expanding perturbatively the Green's function of a system to all orders in 
the impurity scattering. In order to be able to apply the T-matrix approximation one 
needs to have the exact form of the tight-binding Hamiltonian of a system. The T-matrix 
approximation is valid as long as the impurity potential considered is localized, since it is this 
condition that allows the re-summation of all orders in perturbation theory. For extended 
impurities, the T-matrix approximation is in general replaced by the Born approximation, 
for which only the first order term in the impurity potential series is considered; this is 
equivalent to a perturbative expansion in the impurity potential which is valid only when 
the impurity scattering is weak. Nevertheless, both the T-matrix approximation and the 
Born approximation yield the same dependence of the DOS on position, what is different is 
the dependence of the DOS with energy, which is not the main point of our analysis. 

T-matrix approximation in momentum space 



We will briefly review here the principle of the T-matrix approximation |50l . |51| ]. The 
impurity scattering problem can be solved both in the real space and in momentum space. 
We will first focus on the momentum space calculation. For a given system one can define 
a finite temperature (imaginary time) generalized Green's function, 

Gih,k,,T) = -Tre-^(^-^)T,V^fc,(r)V^+(0), (4) 

where K = H — nN, e~^^ = Tre~^^, and T,- is the imaginary time ordering operator. For a 
translationally-invariant, disorder- free system the generalized Green's function defined above 
is non-zero only if ki = k2 (momentum is conserved), and the generalized Green's function 
reduces to the standard Green's function which depends only on one momentum. However, 
the generalized Green's function acquires a non-zero component for ki ^ ^2 if the system 
is inhomogeneous, such as in the presence of impurity scattering. This component can be 



calculated using the T-matrix formulation [58H61|: 

G{ki, k2, iun) = Go{ki, iUn)T{ki, k2, iUn)Go{k2, iun), (5) 

where 

Go{k, iun) = {ii^n - 'Hk)~^, (6) 

is the unperturbed Green's function of the homogenous system, Tik is the Hamiltonian, and 

T{ki,k2,iujn) = V{ki,k2) 

+ J2 Vik,, k')Go{k', tUJn)T{k', k2, lUn). (7) 



As detailed in Fig. HJ this expression stems from a perturbative expansion to all orders 
in the impurity scattering. If only the first term of the expansion is preserved, one recovers 
the first order perturbation theory in the impurity scattering potential, also denoted the 
Born approximation. Oftentimes, one assumes that the impurity scattering potential is 
very close to a delta function so that V is independent of k and k' 

'Himp= [ dxVc\x)c{x)S{x)= f Vc\k)c{k') (8) 

J Jk,k' 

For this case we can solve Eq. ([7]), and obtain 

T(^a;„) = [1-V [ P^Goik, ^co^)]-W (9) 

where Sbz is the area of the first Brillouin zone (BZ) of the system, and the integral over k 
is performed over the entire Brillouin zone. 

In the neighborhood of the impurity, spatial oscillations of the local density of states 
are induced. The Fourier transform of these oscillations can be related to the generalized 
Green's function, 

p{q,(^) = i^9ik,q,u), (10) 

k 

where g{k,q,u) = G{k,k + g, w) — G*{k + q,k,co), and G{k,k + q,u) is the generalized 
retarded Green's function obtained by analytical continuation iun -^ uj + i6 from the Fourier 
transform of the imaginary-time Green's function G{k, k + q, iUn)- 

T-matrix approximation in real space 

If one is interested in the real space spectral properties of a system, such as the space 
dependence of the LDOS, one can use the real space T-matrix formalism. The relations 
described above can be Fourier transformed to real space, yielding for the retarded space- 
dependent generalized Green's function of the system 

G{Ri, R2, E)= f GoiR - Ri, E)T{R, E)Go{R2 - R, E) (11) 

Jii 

where T{R, E) is the space-dependent T-matrix. Same as for the momentum-dependent 
formalism, the generalized Green's function depends on two (spatial) variables. In the 
absence of disorder, for an homogeneous system it is only a function of the difference R1 — R2, 
so that only one spatial index is preserved in the notation: Go{Ri, R2,E) = Go{Ri — R2,E). 
However, in the presence of disorder, the generalized Green's function will depend on both 
variables independently. 

The spatial dependent Green's function characterizes the propagation between two points 
in space Ri and i?2- If one is not interested in propagation of a particle between two spatial 
points, but would rather want to determine the spectral properties, for example the number 
of allowed states having a given energy at a given position, one needs to focus only on the 



limit Ri = R2 = r. Consequently, the LDOS can be obtained from the retarded Green's 
function by using the conventional relation: 

p{f, E) = -lmG(f, f, E) (12) 

In general one focuses on a delta-function impurity localized at i? = 0, which yields 

G(f, f, E) = Go(-f, E)T{E)Go{f, E) (13) 

where T[uj) = [1 — V f -^-^GQ^k^u)]^^]/^, and the integral over k is performed on the first 
BZ, whose area is denoted as Sbz- 

Graphene FT-STS fundamentals 

The Fourier-transform spectroscopy has begun to be widely used also for graphene sys- 
tems. Beyond the fascinating properties of this system and the richness of the possible 
fundamental studies, graphene is for us an interesting system particularly for its band struc- 
ture and for the CEC contour topology at specific energies, which provide a good testing 
ground for the validity of the FT-STS technique. Figure [5] recalls the key feature of the band 
structure of graphene. The CECs (given in a)) start from Dirac points around K points as 
circular contours which become triangular with increasing energy until the contours touch 
together at the M points at the so called Van Hove singularities (VHs). These points are 
indicated in the CEC map, on the 3D view of the band structure (given in b)), and on the 
corresponding density of states. In d) we have schematized the CEC map for an energy near 
the Dirac point and the possible scattering events expected in the presence of impurities. 

Graphene layers have two inequivalent atoms par unit cells. This is one of the key 
properties of the graphene. For this reason graphene should be described by two Bloch 
waves families, and the K points are inequivalent (non-equivalent wave vectors denoted K 
and K'). Around these K points, and close to zero energy, it has been showed that an effective 
Hamiltonian could be defined as 7l|: 



Hk = vpa.p 

where a, denoted pseudospin, is an operator which generates the transformation H^ = 
—Hk'- The pseudospin provides a supplementary quantum number defined only at low en- 
ergy in the immediate vicinity of the K points. The pseudospin is symbolized by the blue 
and red arrows in d). As previously discussed for Erbium disilicide, the expected allowed 
scattering g* wave- vectors could join two points of the same isocontour (intra- valley scatter- 
ing) leading to a circle centered at F, or join two different valleys (inter- valley scattering). 
The first type of scattering leads to long wavelength oscillations in real space, while the 
second one leads to the well known short wavelength reconstruction observed in the vicinity 
of point defects in graphene, and to the \/3 x \/3 reconstruction oftentimes reported for 



graphite and graphene around point defects [72|, |73| and more recently near step edges \M 

However, as we will describe in more details in the next section, not all scattering processes 
are allowed for graphene. Since the pseudospin acts as an extra quantum number, two 
waves of opposite pseudospin cannot interfere and generate standing waves. Thus, since for 
two points on a CEC having opposite momenta relative to a K point, the pseudospins are 
opposite, the intra-valley backscattering is expected to be suppressed. However this is not 
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the case in a backscattering event K, —K connecting two different opposite K valleys. In 
the following we will see how the FT-STS technique could be useful to probe and analyze 
the quantum nature of the quasiparticles on graphene, and how the interpretation of the FT 
can be done using the T-matrix approximation. 

The Hamiltonian and Green's functions for monolayer and bilayer graphene 



The momentum-space tight-binding Hamiltonian for monolayer graphene [75 



is: 



H= (fk[a\h-J{k) + h.c], (14) 

where the operators a\ b^ correspond to creating electrons on the sublattice A and B 
respectively, and f{k) = — t ^ . -|^exp(iA;-aj). Here ai = a{\/3x + y)/2, 02 = a{—\/3x + y)/2, 
03 = —ay, t is the nearest-neighbor hopping amplitude, and a is the spacing between two 
adjacent carbon atoms, which we are setting to 1. 

This form of the Hamiltonian has been used in Ref. [6I1] to perform a numerical analysis 
of the FT-STS spectra. It is also useful to expand the Hamiltonian close to the corners 
of the BZ, which we also denote as nodes or "Dirac points", and use the linearized form 
to solve the problem analytically at low energies [61|. The momenta of the six corners of 
the Brillouin zone are given by i?i,2 = [±47r/(3y3), 0], i?3,4 = [±271/(3^3), 27r/3], K^^q = 
[±27r/(3\/3), — 27r/3]. Close to each corner, m, of the BZ we can write f{q + Km) ~ 4'm{<i) = 
VmQ ■ Jm, where g denotes the distance from the respective corner. Also 111^2 = 3t/2 = v, 
V3^4 = vexp{-in/3), v^^ = vexp{i7i/3) and Ji,2 = (±1, -i), h,i = ^5,6 = (±1,0- 

— * 

The corresponding Green's function, Q(k,uj), derived from the tight-binding Hamiltonian 
in Eq. flT4|) can be expanded at low energy around the six nodes (denoted m), and in the 
2x2 (A,B) sublattice basis can be written as: 

g{k, CO) ^ GUk, uo) = L^ h+^ ~^-(^)\ (15) 

where 6 is the quasiparticle inverse lifetime. The Fourier transform of the linearized Green's 
function is given by (6l|: 



^<pum?{^) h!,'\z) 



G^{f,u;)^u[^.JZ:;,l.^^ ^""^^/ir/ M (16) 



where z = ojr/v, Hq {(r) are Hankel functions, r = |?^, and 0m(^) = Vmr- Jm/{vr). 

On the other hand, the bilayer graphene consists of two graphene layers stacked on top of 
each other such that the atoms in the sublattice A of the first layer occur naturally directly 



on top of the atoms in the sublattice B of the second layer |76l . |77[ , with a tunneling coupling 



of tp. 



In the sublattice basis {A, B) this yields [78|, l79 



Hr-'^(fc) = (^^^J^^^p f'^-f ^'1 (17) 



11 



where for simplicity we have set the effective mass of the quadratic spectrum to 1. The 
corresponding Green's function in real space is given by: 



where we denoted z = ry'\ijj\/v. 

Calculations of the FT-STS for graphene using the T-matrix 



Ref. [61| focuses on monolayer graphene, with a delta-function impurity localized on an 
atom belonging to sublattice A. In the {A, B) basis the impurity potential matrix V has 
only one non-zero component Vn = u. The T-matrix formalism presented above can be 
generalized to graphene for which the Green's functions and the T-matrix for graphene 
become 2x2 matrices, such that 



G{ki, k2, ioJn) = Qo{ki, iUn)T{ki, /cg, i^n)Qo{k2, iun), (19) 



and where 



and T{uj) = [I — V f -^-^Q{k,uj)]~^V, where / is the 2x2 identity matrix, and the integral 

over k is performed on the BZ, whose area is Sbz = Sn'^/Sy/S. 

At arbitrary energy this cannot be calculated analytically, but can be analyzed numeri- 
cally [611 , and the resulting FT-STS spectra (corresponding to the real part of the Fourier 
transform of the LDOS) are plotted in Fig. [61 

The calculation shows regions of high intensity in the FT-STS spectra corresponding to 
intranodal quasiparticle scattering (central region) and internodal scattering (outer regions). 
The central high-intensity region is a filled circle, while the outer regions are empty. Also, 
the rotational symmetry of the high-intensity regions located at the corners of the BZ is 
broken. 

We can compare this result to the one obtained by the joint density of states (JDOS) 
formalism. The JDOS formalism focuses mainly on the position in k-space of the quasipar- 
ticle peaks at a given energy, and considers that scattering takes place equally between all 
quasiparticles living at a given energy. Formally this is equivalent to writing 

p{q,iUn) = lm{TT[Go{k,iuJn)]}Vlm{TT[Go{k + q,iuJn)]}, (21) 

Jk 

We can see that in this formalism neither the chiral structure of the Hamiltonian of 
graphene, nor the phase of the matrix elements are taken into account, and one focuses 
solely on the eigenvalues of the Hamiltonian. The loss of information is evident when in 
Fig. [7] we plot the JDOS result for graphene: 

Indeed here the calculated FT-STS shows a central circle corresponding to the intravalley 
scattering which is not obtained in the T-matrix calculation, nor observed experimentally. 
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The feature observed in the middle of each K contour (inter- valley scattering) is due to the 
trigonal warping of the contour. 



For bilayer graphene, in Ref. [61|, one has considered the case of an impurity located on 
the sublattice A. The resulting FT-STS spectra for the LDOS in the top layer are presented 
in Fig. E 

There are similarities and discrepancies between the monolayer and bilayer cases. Like 
in the monolayer case, there are areas of high intensity centered on the center and corners 
of the BZ. The main difference at low energy is that the central region of high intensity is 
an empty circle, and not a full circle (as for the monolayer case). At high energy, we also 
note a doubling of the number of high intensity lines corresponding to the doubling of the 
number of bands. 



Calculations of the spatial dependence of the LDOS 

We now turn to the study of the dependence of the LDOS on the relative position with 
respect to the impurity (r). We note that a calculation of the LDOS at arbitrary energy 
can be performed numerically using Eqs. (12,13), where the real-space Green's function for 
graphene can be calculated numerically by taking a FT of the full k-space Green's function 
of graphene in Eq. ([H]). However, at low energies, the physics is dominated by linearly 
dispersing quasiparticles close to the Dirac points, and the calculation can be performed 
analytically. The spatial variations of the LDOS due to the impurity have been found to be 
El: 



p{f,E) oc -lm[gi-f,E)T{E)g{f,E)] ^ - 5^Im[e^(^--^")-^G„(-f, E)T(E)G„(f, E)](22) 

m,n 

where m, n denote the corresponding Dirac points. Here T{E) is again the T-matrix, which 
for a delta-function impurity is given by T{u) = [I — V J ■^-^Q{k,u)]~^V, where I is the 
2x2 identity matrix, and the integral over k is performed on the BZ, whose area is Sbz = 

Using Eq. fITB]) and expanding the Hankel functions to leading order in 1/r, in Ref. J6l|, 
it has been found that far from the impurity {ujr/v ^ 1) the corrections to the local density 
of states due to scattering between the nodes m and n are given by: 



p^„(f,a;) oc -Im{t(a;)e^(^--^")-+2W-z[l - 0;(f)0„(r-)] }. (23) 



where t{(jj) is the non-zero element of the T-matrix (Tn), as it was noted that (pi—r) = —(f){r). 

In the case of intranodal scattering {m = n) the above expression vanishes and the LDOS 

is dominated by the next leading correction Pmi'f^-, ^) oc sin{2ur/v)/r'^. This is different from 



what usual wisdom would suggest for a two-dimensional system (1/r decay) [60|, ISOj, and 
has also been described in Refs. 8l|, l82 • 



As briefly outlined above in section Graphene FT-STS fundamentals, the underlying 
physics of this result stems from the chirality of the graphene quasiparticles |83l . |8J] . The 



graphene quasiparticles, due to the presence of two atomic sublattices, have an additional 
degree of freedom deemed pseudospin, such that one says that a quasiparticle belonging to 
sublattice A has pseudospin "up", and a quasiparticle belonging to the sublattice B has 
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pseudospin "down" . By chirality, the pseudospin vector is parallel to the quasiparticle mo- 
mentum, so that if we know the momentum of a quasiparticle, we automatically know its 
pseudospin, thus to which sublattice it belongs and in which proportion. Due to the chiral 
properties of the quasiparticles, backscattering of quasiparticles by extended impurities is 



forbidden [85|, |86|, because when a quasiparticle is backscattered flipping its momentum 
by 180°, by chirality, its pseudospin also flips. Since an extended impurity cannot flip the 
pseudospin of an electron (the wavelengths associated with backscattering are much larger 
than the atomic lattice constant), the backscattering cannot take place. Nevertheless, if 
the impurity is localized, the backscattering is not forbidden 87|. However, the incident 



and backscattered particle have opposite pseudopspins, and thus their wavefunctions can- 
not interfere constructively, same as the wavefunctions of two electrons with spin up and 
spin down cannot give rise to constructive interferences. This lack of interference yields a 
reduction of the Friedel oscillations in the vicinity of the impurity, which no longer decay as 
1/r as in a regular two-dimensional electron gas, but much faster, as 1/r^. 

We should note that the two-dimensional FT of sin(2a;r/f)/r^ is roughly pm{q,(^) oc 
7i9{2u — qv)/2 + arcsin(2c(;/gt')[l — 9{2ijj — qv)]. This corresponds to a filled circle of high 
intensity in the FT-STS spectrum, which is consistent with the results of our numerical 
analysis for the central region of high intensity. 

Nevertheless, for the decay of the Friedel oscillations generated by internodal scattering 
{m ^ n), the chirality considerations are not longer applicable (quasiparticles close to differ- 
ent nodes have different chiralities), and the leading order behavior of the Friedel oscillations 
is 1/r. The FT of cos{2u!r/v)/r is 9{qv — 2ijS) j \/ (fv^ — 4u;^, which translates into empty 
circles of high intensity in the FT-STS spectra, consistent with our numerical analysis. 

For bilayer graphene, an analytical study can be performed at low energies starting from 
the expansion of the Hamiltonian around the Dirac points m. Starting from Eq. (l22l) . 
Ref. [Gli has performed a similar analysis to the case of monolayer graphene. It was noted 
that at large distances {z ^ 1), as opposed to the monolayer case, the leading (1/r) contri- 
bution for intranodal scattering is non-vanishing: 



Pm{r,uj) (X — j=cos{ry/\uj\/v). (24) 

\ui\ 



This is consistent with the appearance of an empty circular contour at the center of the BZ, 
as opposed to the filled circle for the monolayer case. The leading contribution to the decay 
of the oscillations due to internodal scattering is also 1/r. 

Experimental measurements of FT-STS in graphene 

Experiments to measure the LDOS in graphene have been performed by quite a few 
groups J88|492]. Thus, for example, we present below some real-space images from Ref. 90 



for graphene monolayer (Fig. Oa)) and graphene bilayer (Fig. |9lb)). Note that both images 
exhibit a triangular pattern of periodicity ~ 1.9 nm which is related to the interface recon- 



struction [88|, |89|, and which appears as a sextuplet of bright spots in the corresponding 
FFT images([Figs. Oc) andlHld)). Also in Fig. [TUl we present the LDOS and its Fourier 
transform as measured in Ref. |91] . 

The central region of the FFT in Figs. |9lc) and[9ld) is related to intravalley scattering: a 
clear ring-like feature of average radius 1.2 nm~^ is found for the bilayer (Figs. IHd) andlUg)). 
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This radius value is in agreement with the value given in ref. 89|, and with the 2qp value 



derived from ARPES [93|, |9J|. On the monolayer terrace, no central ring is found (Figs. [He) 
and [He)), despite the unprecedented momentum resolution (the result has been checked on 
many different monolayer terraces). This corresponds to the presence of slowly-decaying 
(1/r) long-wavelength oscillations in bilayer graphene, but not in monolayer graphene. 

In the high frequency regions of the FFTs of Figs. [He) and[nid), six outer pockets with 
ring-like shapes centered at K{K') points have been observed. They result from intervalley 
scattering, associated to real space LDOS modulations with a (-\/3 x \/3)R30° periodicity 



with respect to graphene [88|, |89|. As shown in Figs. [Hf) and[9lh), the intensity of the high 
frequency rings in the FFT is not isotropic. The anisotropy is much more pronounced for 
graphene monolayer. The presence of the anisotropy is in agreement with the theoretical 
calculations of Ref. J6ll presented in Fig. [6] 

The experiment described in |90| thus confirms the theoretical picture presented in section 
Calculations of the FT-STS for graphene using the T-matrix, confirming the chiral properties 
of the graphene quasiparticles, and also that these properties can be probed at the nanometer 
scale using scanning tunneling spectroscopy. 



Determination of the band structure of graphene from the FT-STS spectra 

The Fourier transform of the LDOS can be used not only to extract information about 
the chirality of quasiparticles, as described above, but also about the quasiparticle energy 



dispersion (see for example |23l. l89l. l95l. l96l| ) . In Fig.[TT]we present intensity profiles along the 
K — T — K direction for the Fourier transform of the LDOS at various energies |911]. Every 
image was acquired using a lock-in amplifier and a modulation voltage of ±20 meV which 
gave the energy uncertainty. For a wide range of energies, between energies well below the 
Fermi level (—950 meV) to about —150 meV, the width of the peak at the center of the 
Brillouin zone (measured at half maximum) shows a clear dispersion with the bias voltage (it 
decreases when the energy moves towards the Fermi level). At energies closer to the Fermi 
level the profile becomes more complex and displays a dip in the intensity profile close to 
the r point. The central peak also shows two lateral structures (shoulders). The width of 
the central feature does not appear to disperse with energy close to the Dirac point, but it 
starts dispersing again for higher positive energies. 

Fig. [I2] shows a linear dependence of the width of the central feature with energy, except 
in an energy range from —200 meV to 100 meV. For these energies the shape of the central 
peak is more complex. Also, in this range of energy. Fig. [TT]c) presenting the measurement 
of the LDOS as a function of energy taken close to the impurity, shows a 150 meV gap-like 
feature centered around the Fermi level at meV. The Dirac point is estimated at —100 
meV below the Fermi level. The spreading in the dispersion of the points close to the Fermi 
level could be attributed to the presence of the gap. 

The width of the features at the K points follows a similar dispersion. Beyond +500me\^ 



all the features seem to follow a different dispersion branch, as indicated in Fig. [121 In ref. [91 



a comparison between the experimental results with similar theoretical profiles obtained 



IS 



using a T-matrix approximation for a single localized impurity in bilayer graphene [61 
also presented. No gap at the Fermi level was included in the theoretical calculation, while 
a small gap of ^ 100 meV was assumed near the Dirac point, which was taken to be close to 
—250 meV. In Fig. [T21b) is plotted the dispersion of the central and K-points features, and 
of the shoulder, as obtained from the theoretical curves. One observes the presence of two 
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different dispersion branches. The second dispersion branch and the shoulder arise because 
of the bilayer graphene upper band which opens at energies higher than the inter-layer 
coupling. For monolayer graphene these extra features should not appear. 

Thus the measurements in |9l|] show that the quasiparticle approximation and the Fermi 
liquid theory are robust over a large range of energies. The complex structure of the central 
feature (the existence of the shoulder), as well as the presence of two dispersion branches 
are consistent with the bilayer (or multilayer) character of the graphene sample. While the 
point defect modifies the electronic wave-function in its vicinity, a clear linear dispersion 
is still observed, and the relativistic character of the quasiparticle is preserved. Also, the 
STS measurements in Fig. [TTlc) indicate the presence of a gap centered at the Fermi point, 
and not at the Dirac point inferred from our FT-STM measurements. This gap is thus 



different from the Dirac-point gap observed by ARPES ^, |9J] which was attributed to the 
different doping levels of the epitaxial graphene layers. It is however consistent with previous 
STM measurements performed on epitaxial graphene |97j, and more recently on exfoliated 
graphene |98| , where a gap observed at the Fermi level was attributed to a pinning of the 
tunneling spectrum due to the coupling with phonons. 

Measurement of the spatial dependence of the LDOS in the vicinity of a defect 

Figure [131 a) shows a topographic image of a large graphene terrace taken at —17 meV 



(probing filled states), as measured in Ref. [9l[. This layer shows an intriguing "star-like" 



defect with an apparent six-fold {C6v) symmetry. This atomic defect is accompanied by 
a strong distortion of the graphene lattice. The center appears black, which is a dramatic 
change from the case of the unperturbed lattice. As schematized in figure [T3lc), a detailed 
analysis of the real-space image, shows that the point defect is directly located above or 
below a lattice site. 



Using FFT filtering, Ref. [91| has removed the lattice periodicity vectors and all other 
features with wavevectors outside the yellow circle in Fig. [IHlb), thus taking into account 
only the intra-valley and the inter-consecutive-valley scattering processes. The resulting 
real-space image is shown in figure [T51b). This operation strongly enhances the anisotropic 
intensity observed also on the bare topographic image. The LDOS near the defect shows a 
clear threefold (C3f ) symmetry. While, as depicted in Fig. [T3lc), close to the impurity one 
observes a fairly homogeneous standing waves ring, with an almost perfect six-fold symmetry, 
farther away from the impurity, the intensity is clearly higher along three axes (drawn in 
green in Figs. [131 c) and[T3ld)). 

In Ref. [911] the possible origin of this three- fold symmetry has been analyzed by comparing 
the theoretical LDOS for a point defect to the experimental data. The theoretical LDOS 
modulation for a bilayer graphene in the presence of a single impurity has been calculated 
using the real-space T-matrix formalism, modified to incorporate the finite spatial extend 
of the carbon electronic orbitals [99J] . The results for an energy of 300 meV above the Fermi 
level are presented in Fig. [T3ld). One notes quite a few features similar to the ones depicted 
in Fig. [T3lc), including the existence of a three- fold symmetry. However the three-fold 
anisotropy is much less pronounced. 

Figure [m displays the calculated real-space modulations in the LDOS close to the point 
defect, when only specific scattering processes have been considered. For simplicity this 
calculation has been performed for a monolayer graphene, as for bilayer graphene it is less 
clear how much the two bands and the two layers contribute to the measured LDOS. The 
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intra-nodal scattering processes are responsible for radially-symmetric features with small 
q wavevectors and should therefore not be responsible for the presence of the three-fold 
symmetry. Fig. [THh) depicts the sum of the A and B contributions to the full LDOS, 
where in order to increase the anisotropy, the if — )■ i^i B contribution has been multiplied 
by three. This may mimic bilayer graphene where, due to the coupling between the two 
layers, the sub lattices A and B may not contribute equally to the observed LDOS. A very 
strong threefold symmetry is observed in this weighted superposition, and the resulting 
image corresponds more closely to the pattern observed experimentally in figure [T3lb). This 
is consistent with having a defect at an A site whose dominant effect is the scattering of 
sublattice B electrons between two conse cutiv e valleys [K — )■ K\). 

More recent experiments presented in [lOO] indicate that in monolayer graphene patterns 
similar to the ones observed in [91] but six-fold symmetry may arise as a result of rotating 
sequences of dislocations that close on themselves, forming grain boundary loops that either 
conserve the number of atoms in the hexagonal lattice or accommodate vacancy/interstitial 
reconstruction, whil e lea ving no unsatisfied bonds. It would be interesting to generalize the 



results presented in llOOl ] to bilayer graphene to see if one can obtain similar patterns to the 
ones observed in 19 ll. 



Van Hove extension deduced from FT-STS features on graphene modified by the 
intercalation of gold clusters. 

The FT-STS technique has recently also been used to interpret the strong standing waves 
pattern observed on epitaxial monolayer (ML) graphene, modified by the intercalation of 
gold atoms. As described also above, epitaxial graphene on SiC(OOOl) consists of a buffer 
graphene layer which is covalently bonded with the substrate, and of a ML graphene weakly 
connected to the buffer layer. The monolayer graphene is n-doped with a transfer of electrons 
from the substrate. We have discovered that the deposition of gold atoms under UHV at 
room temperature, followed by an annealing cycle, leads to the intercalation of gold in 



different forms [36[. One of them is the intercalation of aggregates made of 1 to 3 small flat 
6-atoms gold clusters. Figure [15] shows these results. In a) is shown a topographic image 
taken at -IV. The clusters appear under the graphene ML as bright spots. The dl/dV map 
image of the same area taken at +0.9V shows that the clusters appear as dark regions. As 
we probe here the empty states it seems that the aggregates create an excess of electron as 
schematized in d). 

In between the gold clusters, a standing waves pattern develops, as revealed by bright 
p(2x2) protrusions, and the size of the dark regions decreases as the bias voltage increases 



92[. These standing waves are associated to elliptic features in the FT-STS spectra, located 
around the M points as shown in c). The size of these elliptic features increases linearly 
with the bias voltage. Using the JDOS calculation we have deduced that these features are 
associated to the VHs (92|. In fact these elliptic features are observed in the JDOS spectra 
only when the CEC touch each other at the M point. This is shown in the calculated JDOS 
e) and h) for the CEC f ) and g) respectively. In f ) the contours do not touch each other, 
and no features are observed around M in e) while for h) as soon as the contours touch each 
other, an ellipse is observed. 

This ellipse is associated to the k-vectors symbolized as blue arrows, which connect the 
apex of two consecutive triangular K contours (inter- valley scattering). The size of the 
ellipse is associated to the filling of the states located at the triangle apex. This indicates 
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an extension of the V Hs. T his has been very recently confirmed by ARPES measurements 
in the occupied states [l01|. Two other features in the perpendicular direction T — M — T 
expected in both JDOS, and even in T-matrix approximation calculations are not observed 
experimentally. These features are associated to the wavevectors symbolized by the red 
arrows which should correspond to intra-valley scattering. As we are at large energies 
this could not be due to the pseudospin orientation. The possibility to a nodal- ant inodal 
dichotomy has been discussed but remains under debate [92!]. The reasons why the clusters 
create such an VHs extension remains non understood. We have recently shown that a 
change in the third nearest neighbor hopping energy in the tight-binding Hamiltonian create 
strong modifications of the band structure around the M points and even give rise to new 



Dirac points |102| . but the experimentally observed third- nearest neighbor coupling seems 
at present too small to justify important modifications of the band structure for our system. 
More experimental studies and ab-initio calculations are currently undergoing to clarify the 
nature of these standing waves. 



GENERAL CONCLUSION 

We have shown that simple JDOS calculations and their comparison with the FT-STS 
maps provide an accurate determination of the size and shape of circular, and even non- 
circular free electron-like constant energy contours. These contours in turn provide precise 
information about the quasiparticle dispersion, and fine details about fairly complex 2D band 
structures. On the other hand, we have shown that the T-matrix approximation becomes 
absolutely necessary for more complex systems where more than one quantum number is 
involved in the scattering process. 

While in order to apply the T-matrix approximation one needs to have the exact form of 
the tight-binding Hamiltonian of a system, the JDOS approach does not rely on the exact 
form of the Hamiltonian, but on the phenomenological form of the equal energy contours. 
Its advantage is that it can be used even if the Hamiltonian of a system is not known, 
but the equal energy contours are known from a different experiment such as ARPES; this 
allows one to have an intuitive reading of the FT features. The results of the T-matrix and 
the JDOS approximations are in general the same for a symmetrical, quadratic-dispersing 
system; for example, for the case of Au(lll), while not shown here, both the T-matrix and 
the JDOS approximations yield circles of high intensity in the FT-STS spectra. However, 
for more complicated systems, these two methods may yield different results, and one notes 
that when the features predicted by JDOS are not observed experimentally, as it is the case 
of graphene, the physics of the system is in general more complicated, and a full T-matrix 
calculation is needed to describe it. 

We notice that in general a good agreement with angle resolved photoemission data has 
been found and this indicates that both methods (ARPES and FT-STS) probe the same 
quasiparticle excitations of the system. This technique provides the possibility to measure 
the band dispersion modification generated by specific defects of different nature. This 
allows to test locally the effect of specific modifications and functionalization of well-known 
surfaces which is not possible with ARPES measurements which requires to prepare large 
homogeneous surface in order to obtain reliable data. 
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Figures 
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FIG. 1: The STM picture a) shows the herringbone reconstruction of the Au(lll) surface island, 
perturbed by a nearby dislocation (50 x 50nm, -0.4 V). b) is the dl/dV map of the topographic 
image a). Three 2D-FFTs of these dl/dV maps (with the corresponding applied bias) are shown 
in c), from which we get the energy dispersion curve of Au(lll) shown in d). f) Topographic and 
g) dl/dV map at -30 meV (20 x 20nm) showing both standing waves and atomic resolution on 
ErSi2 islands on Si(lll). Three 2D-FFTs of these dl/dV maps (with the corresponding applied 
bias) are shown in h), from which we get the energy dispersion curve of ErSi2/Si{lll) shown in 
i). e) and j) represent the geometrical construction to get access to the density of scattering events 
of a CEC for gold and erbium disilicide respectively. 
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FIG. 2: Constant Energy Contours of ErSi2 in the energy region where both hole and electron 
band coexist (-200 meV to + 200 nieV). b) shows the experimental band structure of ErSi2 
determined with ARPES and FT-STS measurements. The band structure is compared to DFT 
calculation of YSi2- c) from left to right : Fourier transform of an STM image of ErSi2- The 
dominating structures observed in M can be modelized by JDOS construction as depicted in Fig[TJ 
JDOS calculations showing the contribution of different scattering events: the third frame shows 
the scattering events between the hole-like band and the electron-like band whereas the fourth 
frame displays the structure obtained by only electron-like band scattering. 
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FIG. 3: a) Representation of the stationary phase (blue parallel vectors) of a scattering event (red 
vector) resulting from the scattering from one state to another (depicted by the grey wave vectors) 
42l |. b) CEC from erbium disilicide, and in c) the corresponding JDOS calculation, d) shows 
the dominating structures from c), when the stationary phase condition is taken into account; the 
result is in good agreement with the experimental FFT from Fig. [21a). 
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FIG. 4: Diagrammatic description of the T-matrix approximation: the Green's function for a 
system in the presence of an impurity can be obtained from adding to the unperturbed Green's 
function the contributions corresponding to all-order impurity scattering processes. We note that 
the sum of the contributions appearing to the right of Go is identical to T, which allows one to 
write the self-consistency equation on the third line. 
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FIG. 5: A three-dimensional (3D) representation of the band structure of graphene, and a plot 
of the constant-energy contours, a) The constant energy contours starting from the Dirac point, 
calculated in the nearest-neighbor tight-binding approximation. The key features of the band 
structure, the Dirac points and the Van Hove singularities are indicated on the 3D representation 
in b), and on the theoretical density of state in c). The expected intra- and inter- valley scattering 
processes and the corresponding pseudospin orientations (as discussed in the text) are schematized 
in d) 
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a) 



b) 



FIG. 6: FT-STS spectra for a monolayer graphene sample with a single delta-function impurity. 
Fig. a) corresponds to energies 0.15t at 5 = 0.07t. The actual lowest (0) and highest (1) values 
of the FT-STS intensity are (—1.3,2.6) in arbitrary units. Fig. b) shows a cross section of the 
FT-STS intensity as a function of kx for ky = for energy 0.15t. 




FIG. 7: In b) the FT-STS map for a monolayer graphene calculated using the JDOS approximation 
for the CEC given in a). The arrows in b) show the inter- valley scattering momenta leading to 
the feature observed around the K points in b). The central circle corresponding to intra- valley 
scattering does not appear in the T-matrix calculation, and is not observed experimentally. As 
detailed in the text, this is a consequence of the chirality of graphene quasiparticles. 
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FIG. 8: FT-STS maps for a bilayer sample. Figs, a) and c) depict the FT-STS intensity in arbitrary 
units at energies O.li, 0.4t, with tp = 0.3t, and 6 = 0.05t. The actual lowest (0) and highest (1) 
values of the FT-STS intensity are (—6.9, 3.6) and (—6.0, 7.4) respectively. Figures b) and d) are 
cross-sections at ky = of Figs, a) and c). 
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FIG. 9: a),b) Low-bias STM images of 50 nm wide monolayer a) and bilayer b) terraces. Sample 
bias and tunneling current are respectively +2 mV and 0.4 nA for a), +4 mV and 0.13 nA for 
b). c), d) Two-dimensional fast Fourier transform (FFT) maps of the STM images a) and b). e) 
Central region of c), showing no intravalley-backscattering related ring (the green arrow points out 
the position where such a ring should appear), f) One of the outer pockets of c). g) Central region 
of d), showing a clear ring- like feature of radius 2qF related to intravalley-backscattering. h) One 
of the outer pockets of d). Outer pockets shown in f) and h) are centered at the K (or K^) point 
and result from intervalley scattering. Data courtesy : I. Brihuega et al. |9C 





11 




'"'i*^*« 




•••'••• 


'•■'■■'■' 




FIG. 10: a) STM topographic image (10 x lOnm^, -17meV, InA) showing an isolated defect, b) 
The FFT power spectrum of the 2D topographic image in a). The features indicated by the arrows 
correspond to inter-valley scattering. 
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FIG. 11: a), b) Intensity profiles along the K — T — K direction of tlie 2D FT of the LDOS taken 
for different bias voltages, c) STS measurement taken near the impurity showing a gap centered 
around the Fermi level. 
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FIG. 12: The dispersion for the width of the central-ring feature (red circles) and of the K-points 
feature (black triangles) obtained in a) from the experimental data from Fig. \TT\ and from the 
theoretical K — T — K cuts (b)). The dispersion of the central- feature shoulder on the theoretical 
curves is indicated by blue squares. 
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FIG. 13: a) STM topographic image (10 x lOnm^, -UmeV, InA) showing an isolated defect, b) 
Zoom- in (5 x ^nm?) on the defect, subsequent to a FFT filtering removing atomic resolution and 
second order features in the 2D FFT. c) 2.5 x 2.5nm'^ zoom in on b) with the schematic atomic 
lattice of graphene (A/black and B/yellow atoms) superimposed over the standing waves pattern 
around the point defect, d) The real-space LDOS in a bilayer graphene calculated for a point 
defect using a single-impurity T-matrix approximation at SOOmeV above the Dirac point. 
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FIG. 14: Calculated real-space LDOS modulations in a monolayer graphene for a point defect 
placed on top of an A atom. Figs. 14 a), b), c) depict the contributions of selected scattering 
processes to the LDOS {K — )■ Ki,K — )■ K2 and K — )■ K^ respectively) as indicated in Fig. 5 d); 
Figs. 14 f), and g) depict the separate contributions of the A and B sublattices to the K — )■ Ki 
LDOS modulations in a). Figs. 14 d) and e) depict the A and B sublattice contributions to the full 
LDOS, evaluated when all scattering processes are considered. In h) we depict the sum of the A 
and B contributions to the full LDOS, when the weight of the K — )■ Ki B component is increased 
three times compared to that of the K — )■ Ki A component. 
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FIG. 15: a) Topographic image taken at -1.4V (13 x 13nm^) of ML graphene with the intercalated 
gold clusters underneath, b) Same region dl/dV map at +0.9V (unoccupied states), c) FT of 
b). e) and h) are the JDOS calculation of the CEC given in f) and g) respectively, i) is a zoom 
of the ellipitic feature encircled in h). Blue and red arrows correspond to the momentum vectors 
associated to the black vectors. The expected features in the JDOS calculation that are not 
observed experimentally are associated to momentum symbolized by red arrows. 
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